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Abstract 

We  establish  rigorously  the  fundamental  nesting  behavior  of  T-spline  spaces  in  terms  of 
the  topology  of  the  T-mesh.  This  provides  a  theoretical  foundation  for  local  refinement 
algorithms  based  on  analysis-suitable  T-splines  and  their  use  in  isogeometric  analysis.  A 
key  result  is  a  dimension  formula  for  smooth  polynomial  spline  spaces  defined  over  the 
Bezier  mesh  of  a  T-spline. 

Keywords:  T-splines,  isogeometric  analysis,  local  refinement,  T-spline  spaces 


1.  Introduction 

T-splines  were  originally  introduced  as  a  superior  alternative  to  NURBS  [1]  and  have 
emerged  as  an  important  technology  across  several  disciplines  including  industrial,  ar¬ 
chitectural,  and  engineering  design,  manufacturing,  and  engineering  analysis.  T-splines 
can  model  complicated  designs  as  a  single,  watertight  geometry  and  can  be  locally  re¬ 
fined  [2,  3].  These  basic  properties  make  it  possible  to  merge  multiple  NURBS  patches 
into  a  single  T-spline  [4,  1]  and  any  trimmed  NURBS  model  can  be  represented  as  a 
watertight  T-spline  [5]. 

T-splines  are  an  ideal  discretization  technology  for  isogeometric  analysis  [6,  7,  3,  8,  9, 
10,  11].  Isogeometric  analysis  was  introduced  in  [12]  and  described  in  detail  in  [13].  The 
isogeometric  paradigm  is  simple:  use  the  smooth  geometric  basis  as  the  basis  for  analysis. 
Traditional  design-through-analysis  procedures  such  as  geometry  clean-up,  defeaturing, 
and  mesh  generation  are  avoided.  Additionally,  the  higher-order  smoothness  provides 
substantial  gains  to  analysis  in  terms  of  accuracy  and  robustness  of  finite  element  solu¬ 
tions  [14,  15]. 

Analysis-suitable  T-splines,  a  mildly  restricted  subset  of  T-splines,  are  optimized  to 
meet  the  needs  of  both  design  and  analysis  [16,  3].  Analysis-suitable  T-splines  main¬ 
tain  the  important  mathematical  properties  of  NURBS,  such  as  linear  independence  [16] 
and  partition  of  unity  [17],  while  providing  an  efficient  and  highly  localized  refinement 
capability  [3]. 
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In  this  paper  we  establish  rigorously  the  nesting  behavior  of  analysis-suitable  T- 
splines.  This  provides  the  theoretical  justification  for  the  analysis-suitable  local  refine¬ 
ment  algorithm  in  [3]  and  provides  important  insights  into  T-spline  spaces  in  general. 
A  key  result  is  a  simple  dimension  formula  for  polynomial  spline  spaces  defined  over 
the  Bezier  mesh  of  a  T-spline  [18].  The  formula,  written  only  in  terms  of  topological 
quantities,  “bridges  the  gap”  between  T-splines  and  traditional  spline  spaces  posed  over 
partitions  of  R2 .  We  feel  that  this  may  have  important  implications  in  establishing  ap¬ 
proximation,  stability,  and  error  estimates  for  T-splines  in  isogeometric  analysis  [19]  by 
allowing  theoretical  tools,  established  for  traditional  spline  spaces,  to  be  used  directly 
for  T-spline  spaces. 

This  paper  is  organized  as  follows.  Section  2  establishes  notation  and  reviews  fun¬ 
damental  T-spline  concepts.  T-junction  extensions  and  the  extended  T-mesh  are  then 
described  in  Section  3.  Section  6  uses  the  smoothing  cofactor-conformality  method  to 
develop  an  algebraic  representation  of  the  smoothness  properties  of  a  polynomial  spline 
space  posed  over  the  Bezier  mesh  of  a  T-spline.  A  dimension  formula  for  these  polyno¬ 
mial  spline  spaces  is  then  derived  in  Section  7.  Finally,  Section  8  establishes  the  nesting 
behavior  of  T-spline  spaces. 

2.  T-spline  fundamentals 

We  restrict  our  developments  to  cubic  T-spline  surfaces.  All  results  can  be  easily 
generalized  to  other  odd  polynomial  orders. 

2.1.  The  T-mesh 

A  T-spline  is  defined  in  terms  of  a  control  grid  or  T-mesh,  T,  and  two  global  knot  vec¬ 
tors,  S  =  [£-2,£-i,£o,---,£c+i,£c+2,£c+3]  and  II  =  [r;_2,  V-i,  Vd+i,  Vd+2,  Vd+sil- 

Interior  knots  may  have  a  multiplicity  of  three  while  end  knots  may  have  a  multi¬ 
plicity  of  four.  The  global  knot  vectors  define  a  full  parametric  domain,  H  C  M2, 
where  fl  =  [£-2,63+3]  <8>  [v-2,ild+3]  and  a  reduced  parametric  domain,  Cl  C  Cl,  where 
Cl  =  [£i,£c]  ®  [Vi,Vd\- 

Each  control  point  or  active  vertex ,  t!“?,  corresponds  to  a  unique  pair  of  knots,  (£, ,  rjj), 
0<i<c+l,  0<j<d+l.  Note  that  active  vertices  may  be  positioned  at  some,  but 
not  all,  pairs  of  indices  (i,j)  within  0  <  i  <  c  +  1,  0  <  j  <  d  +  1  because  of  T-junctions. 
A  T-junction  terminates  a  row  or  column  of  vertices  and  edges  before  the  boundary  of 
the  T-mesh  is  reached.  The  vertices  Vij  =  (£i,rjj),  —2<i<—lorc+2<i<c+3, 
and  —2  <  j  <  —1  or  d  +  2  <  j  <  d  +  3  are  called  inactive  vertices  and  are  not  associated 
with  T-spline  blending  functions.  We  denote  the  total  number  of  vertices  in  T  by  n  and 
the  number  of  active  vertices  by  na .  In  the  following  we  also  specify  a  global  index  A  for 
each  active  T-mesh  vertex  such  that  A  =  r(i,j),  where  r(i,j)  is  an  index  map. 

Figure  1  shows  a  T-mesh.  The  global  indexing  corresponding  to  the  global  knot 
vectors  is  shown  along  the  bottom  and  left.  The  active  vertices  are  denoted  by  open  and 
red  circles  where  the  red  circles  are  T-junctions.  Shaded  circles  denote  inactive  vertices. 
The  full  parametric  domain  is  the  union  of  the  light  and  dark  shaded  rectangles.  The 
reduced  parametric  domain  is  denoted  by  the  dark  shaded  rectangle. 

Each  vertical  edge  is  associated  with  an  6,  —2  <i<c  +  3,  and  each  horizontal  edge 
is  associated  with  an  rjj,  —2  <  j  <  d  +  3.  Two  perpendicular  edges  can  only  intersect 


2 


d+3 

d+2 

d+1 

d 


3 
2 

1 

0 
-1 
-2 

-2-10123  c  c+1 c+2c+3 

Figure  1:  A  T-mesh.  The  global  indexing  is  shown  on  the  bottom  and  left.  The  active  vertices  are 
denoted  by  open  and  red  circles.  Inactive  vertices  are  denoted  by  shaded  circles.  T-junctions  are 
denoted  by  red  circles.  The  full  parametric  domain  is  denoted  by  the  union  of  the  light  and  dark  shaded 
rectangles.  The  reduced  parametric  domain  is  denoted  by  the  dark  shaded  rectangle. 


at  a  vertex.  Each  edge  is  assigned  a  knot  interval  which  is  the  parametric  length  of  the 
edge.  An  edge  segment  is  a  row  or  column  of  T-mesh  vertices  and  edges  which  begins  and 
ends  at  a  T-junction  or  the  boundary  of  the  T-mesh.  We  denote  a  vertical  edge  segment 
by  e\  and  a  horizontal  edge  segment  by  e^.  If  the  orientation  of  the  edge  segment  is 
not  important  we  simply  write  e*,  without  a  superscript.  We  denote  the  number  of  edge 
segments  in  T  by  nseg ,  the  number  of  horizontal  edge  segments  by  nhseg ,  and  the  number 
of  vertical  edge  segments  by  nvseg . 

We  call  each  polygonal  face  in  a  T-mesh  a  T-mesh  element.  The  parametric  domain 
of  each  T-mesh  element  is  denoted  by  tte  C  where  e  is  a  global  element  index.  Note 
that  not  all  T-mesh  elements  are  contained  in  the  reduced  parametric  domain. 

The  notation  T1  C  T2  will  indicate  that  T2  can  be  created  by  adding  vertices  and 
edges  to  T1  and  the  notation  T1  C  T2  will  indicate  that  T2  can  be  created  by  adding  edges 
and  vertices  to  existing  edge  segments  in  T1  without  creating  additional  edge  segments. 
Obviously,  T1  C  T2  implies  T1  CT2. 

2.2.  T-spline  spaces 

A  T-spline  blending  function,  NA(£,  if),  is  associated  with  each  active  T-mesh  vertex. 
The  T-spline  blending  functions  are  given  by 

NA{^r})  =  B[~AmB\nA]{rj)  (1) 
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where  B[Ea]{£.)  and  i?[II^4] (77)  are  the  cubic  B-spline  basis  functions  associated  with  the 
local  knot  vectors 


-A  =  [£;(a,i),  •  •  •  ,&(A,5)]  C  5  (2) 

11a  =  [vj(A,i),---,Vj(A,5)\  cn  (3) 

and  i{A,i)  and  j(A,j )  map  the  local  knot  indices,  i  and  j,  to  the  global  knot  indices, 
i  and  j.  The  knot  vectors  2a  and  IIa  are  inferred  from  the  T-mesh  according  to  the 
procedure  outlined  in  [1],  Figure  2  illustrates  the  construction  of  a  T-spline  blending 
function  corresponding  to  active  vertex  v33.  In  this  case,  the  local  knot  vectors  are 
SA  =  [£1,6,6,  £4,63]  and  ru  =  [171,172, 173,174,175] ■ 
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Figure  2:  Inferring  a  T-spline  blending  function  from  a  T-mesh. 


The  set  of  all  T-splines  with  the  same  T-mesh  topology,  T,  and  global  knot  vectors  is 
called  a  T-spline  space.  We  denote  a  T-spline  space  by  T.  While  the  T-spline  blending 
functions  are  defined  using  the  full  parametric  domain,  12,  a  T-spline  space  is  only  defined 
over  the  reduced  parametric  domain,  12.  In  other  words, 

T  =  {/  S  L2{p.)  \f  =  J2  caNaW,  ca  €  (4) 

where  L2(I2)  is  the  space  of  square  integrable  functions  over  12. 

3.  Extended  T-spline  fundamentals 

To  develop  the  needed  mathematical  machinery  to  describe  a  T-spline  space,  T,  we 
define  the  extended  T-mesh,  T ext,  and  elemental  T-mesh,  T eiem.  Extended  and  elemental 
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T-meshes  are  defined  using  T-junction  extensions. 

3.1.  T-junction  extensions 

A  T-junction  extension  is  a  closed  line  segment  which  begins  at  a  T-junction  and 
extends  into  the  T-mesh  away  from  the  T-junction.  There  are  two  types  of  T-junction 
extensions:  face  and  edge  extensions.  An  n-bay  face  extension  is  a  line  segment  which 
originates  at  a  T-junction  and  extends  in  the  direction  of  a  missing  edge  until  n  perpendic¬ 
ular  edges  or  vertices  are  intersected.  T-junction  face  extensions  include  the  intersection 
points  but  not  the  originating  T-junction.  For  simplicity,  the  term  face  extension  will 
denote  a  2-bay  face  extension.  An  edge  extension  can  be  created  if  an  edge  is  attached  to 
the  T-junction  in  the  opposite  direction  of  a  face  extension.  The  edge  extension  begins 
at  the  T-junction  and  extends  to  the  edge’s  opposite  vertex.  Since  T-junction  extensions 
(face  or  edge)  are  closed  line  segments,  a  horizontal  and  vertical  extension  can  intersect 
either  on  the  interior  of  both  extensions  or  at  the  endpoint  of  one  extension  or  both 
extensions.  Figure  3  shows  several  examples  of  intersecting  T-junction  extensions.  The 
edge  extensions  are  the  dashed  red  lines  and  the  face  extensions  are  the  dotted  black 
lines. 


Figure  3:  Examples  of  T-junction  extension  intersections.  The  edge  extensions  are  the  dashed  red  lines 
and  the  face  extensions  are  the  dotted  black  lines. 


3.2.  The  extended  T-mesh 

An  extended  T-mesh,  Text,  is  a  T-mesh  which  is  formed  by  adding  face  (and  possi¬ 
bly  edge)  extensions  to  T.  We  denote  the  number  of  T-mesh  vertices  in  Text  by  next- 
The  vertices  corresponding  to  crossing  face  extensions  are  called  crossing  vertices.  The 
vertices  corresponding  to  overlapping  face  extensions  are  called  overlap  vertices.  The 
inactive  vertices  and  the  vertices  corresponding  to  the  intersections  of  face  extensions 
and  existing  edges  in  T,  which  do  not  already  correspond  to  overlap  vertices,  are  called 
extended  vertices.  We  denote  the  number  of  crossing,  overlap,  and  extended  vertices  by 
n+ ,  n~ ,  and  n* ,  respectively.  We  note  that  next  =  na  +  n+  +  n~  +  n*  where  na  is  the 
number  of  active  vertices  in  T. 

Figure  4  shows  a  T-mesh  and  the  corresponding  extended  T-mesh,  Text.  The  face 
extension  edges  which  have  been  added  to  T  are  denoted  by  dotted  black  lines.  The  cross¬ 
ing  vertices  are  denoted  by  red  stars,  the  overlap  vertices  are  denoted  by  red  hexagons, 
and  the  extended  vertices  are  denoted  by  red  squares.  Notice  that  the  inactive  vertices 
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are  also  extended  vertices.  The  active  vertices  are  denoted  by  hollow  circles  and  the 
T-junctions  are  denoted  by  red  circles. 

Every  edge  segment  in  Text  is  composed  of  at  least  four  edges  and  five  vertices.  All 
vertices  in  an  edge  segment  which  are  not  extended  vertices  are  called  interior  vertices. 
We  say  an  edge  segment  ei  is  prior  to  edge  segment  e 2  if  there  exists  an  extended  vertex 
in  e2  which  is  an  interior  vertex  in  e\.  A  set  of  edge  segments  *  =  1,2, . . . ,  c  form  a 
circle  if  e,  is  prior  to  e^+i  for  i  =  1, ...  ,c  —  1  and  ec  is  prior  to  ei.  The  intersection 
vertices,  14,  of  the  edge  segments  e^-i  and  e*  in  a  circle  are  called  circle  vertices.  Notice 
that  every  edge  segment  in  a  circle  must  contain  at  least  two  active  vertices. 

3.3.  The  elemental  T-mesh 

T-mesh  edges  may  not  correspond  to  all  knot  lines  in  the  underlying  T-spline  blending 
functions  due  to  T-junctions.  The  elemental  T-mesh,  T e;em,  is  an  extended  T-mesh 
formed  by  adding  the  missing  edges  and  vertices  to  T  such  that  all  blending  function  knot 
lines  are  represented.  The  missing  edges  are  added  in  the  form  of  n-bay  face  extensions 
where  n  is  determined  by  the  knot  structure  of  the  underlying  blending  functions.  All 
notation  and  properties  described  in  Section  3.2  for  the  extended  T-mesh  apply  to  an 
elemental  T-mesh.  Note  that  T  C  Text  and  Text  U  T eiem. 

4.  Perturbing  T-spline  spaces 

It  will  often  be  necessary  to  perturb  a  T-spline  space.  To  perturb  a  T-spline  space 
a  set  of  knot  intervals  is  modified  by  a  small  parameter.  In  a  smoothing  perturbation 
all  zero  knot  intervals  in  a  T-mesh  are  replaced  by  an  arbitrarily  small  knot  interval, 

5.  We  denote  a  smoothed  T-mesh  by  T(<5),  a  smoothed  extended  T-mesh  by  Text(S),  a 
smoothed  elemental  T-mesh  by  Tejem(<5),  and  a  smoothed  T-spline  space  by  T(6).  Notice 
that  by  replacing  all  zero  knot  intervals  with  a  small  parameter  the  resulting  space  is 
globally  C2-continuous.  If  a  smoothing  perturbation  is  performed  on  a  T-mesh  with  no 
zero  knot  intervals  the  T-mesh  is  left  unchanged. 

In  an  offset  perturbation  all  face  extensions  in  Text,  which  overlap,  are  offset  by  an 
arbitrary  small  real  parameter,  e.  We  denote  an  e-offset  T-mesh  by  T(e),  an  e-offset 
extended  T-mesh  by  Text(e),  an  e-offset  elemental  T-mesh  by  Te;em(e),  and  an  e-offset 
T-spline  space  by  T(e).  Notice  that  an  e-offset  extended  T-mesh  does  not  have  any 
overlap  vertices.  If  an  offset  perturbation  is  performed  on  a  T-mesh  which  does  not  have 
any  overlap  vertices  the  T-mesh  is  left  unchanged. 

Since  T-spline  blending  functions  are  continuous  functions  of  their  defining  knot  in¬ 
tervals  we  have  that  T(S)  — >•  T  as  S  — >  0  and  T(e)  — >  T  as  e  — >•  0. 

5.  The  extended  and  elemental  spline  spaces 

Using  the  smoothed  extended  T-mesh,  Text(5),  we  define  the  homogeneous  extended 
spline  space  as 

Sext  =  {/  6  C2’2(R2)  I  e  P33,Vf^t  C  fi,  and/|R2^  =  o}  (5) 

where  U2’2(R2)  is  the  space  of  bivariate  functions  which  are  C2-continuous  in  £  and 
77  over  all  of  R2.  P33  is  the  space  of  bicubic  polynomials.  The  homogeneous  elemental 
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Figure  4:  A  T-mesh  and  the  corresponding  extended  T-mesh.  The  face  extension  edges  which  have  been 
added  to  T  are  denoted  by  the  dotted  black  lines.  The  crossing  vertices  are  denoted  by  red  stars,  the 
overlap  vertices  are  denoted  by  red  hexagons,  and  the  extended  vertices  are  denoted  by  red  squares. 
Notice  that  the  inactive  vertices  are  also  extended 'jvert ices.  The  active  vertices  are  denoted  by  hollow 
circles  and  the  T-junctions  are  denoted  by  red  circles. 


spline  space,  Seiem,  can  be  defined  similarly  using  Te;em(<5).  Notice  that  the  homogeneous 
spline  spaces  are  defined  over  all  of  R2.  Restricting  the  domain  of  definition  generates 
the  extended  spline  space,  Text,  where 

Text  =  (6) 

The  elemental  spline  space,  Tiem,  can  be  defined  similarly  using  Seiem. 

Theorem  5.1.  For  any  T-mesh,  T,  we  have  that 

dim  Text  =  dim  Sext  (7) 

and  similarly  for  Seiem  and  Teiem  ■ 

Proof.  We  first  prove  that  the  dimension  of  Sext  is  not  less  than  the  dimension  of  Text- 
Notice  that  for  any  function  /  £  Sext,  /|q  S  Text-  We  now  show  that  the  dimension  of 
Text  is  not  less  than  the  dimension  of  Sext-  This  is  equivalent  to  showing  that  there  is 
only  one  function  in  Sext  which  is  zero  over  fl.  It  is  easy  to  see  that  the  only  function 

which  is  zero  over  Q  must  be  zero  over  all  of  P2  since  the  minimum  support  of  a  cubic 

C 2  spline  function  is  four  intervals.  □ 

Remark  5.2.  While  we  develop  our  theoretical  results  using  Text(<5),  Sext,  and  Text  the 
results  also  hold  */Te/em(<5),  Seiem,  and  Teiem  are  used  instead. 

6.  The  smoothing  cofactor-conformality  method 

The  smoothing  cofactor-conformality  method  [20]  can  be  used  to  transform  the  com¬ 
plicated  smoothness  properties  of  Sext  into  a  linear  constraint  matrix,  M.  This  constraint 
matrix  can  then  be  analyzed  to  determine  the  dimension  of  Sext-  For  additional  appli¬ 
cations  of  this  method  in  the  context  of  spline  spaces  over  T-meshes  see  [18,  21,  22,  23, 
24,  25,  26]  and  references  therein. 

6.1.  Vertex  and  edge  co factors 

As  shown  in  Figure  5,  for  any  vertex,  Vij  =  ( £,i,r]j ),  in  T ext(5),  the  surrounding 
bi-cubic  polynomial  patches  are  labeled,  pT (£,  r/) .  k  =  0, 1,2,3.  If  the  vertex,  Vij,  is 
a  T-junction,  then  p*  •(£,  ? 7)  =  for  some  k.  Since  Pij(^,v)  and  Pi.j(€ ,v)  are 

C2-continuous  there  exists  a  cubic  polynomial  \T(r]),  called  the  edge  cofactor,  such  that 


Pi,j (£,  v)  ~  Pij (£,  p)  =  A2  j  (7?) (£  -  £i)3 ■  (8) 

Similarly,  there  exists  cubic  polynomials,  Xjj(rj),  p,V(£),  and  pfj(f),  such  that 

Pij  ($,V)~  Pi,j  (£,»?)  =  Pi,j  (9) 

V)  -  Pi,j(t,  v)  =  ~  &)3,  (10) 

p°i,j(£,v) -Pi,j(£,v)  =  ~  Vj)3-  (11) 
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Figure  5:  The  smoothing  cofactors  around  a  vertex. 


We  note  that  if  two  patches  are  identical  the  edge  cofactor  is  zero.  Combining  (8)  -  (11) 
gives 

(AL0?)  -  -  &)3  =  -  mf-  (12) 

Since  (£  —  £;)3  and  (7?  —  rjj)3  are  prime  to  each  other  there  exists  a  constant,  dij,  called 
the  vertex  cofactor ,  such  that 

A,j(v)  ~  xh(v)  =  di,j{r)  -  Vj)3,  A,j( 0  -  M?d(0  =  -  &)3-  (13) 

d.£.  Assembling  the  constraint  matrix,  M 

A  global  constraint  matrix,  M,  can  be  constructed  by  examining  the  coupling  between 
edge  cofactors  on  all  horizontal  and  vertical  edge  segments  in  Text{$). 
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Figure  6:  The  smoothing  cofactors  along  a  horizonal  edge  segment. 


Referring  to  Figure  6,  consider  a  horizontal  edge  segment  with  k  +  1  vertices  and  k 
edge  cofactors.  Using  (13)  we  have  that 

Hioj-o  =  <At-Sio)3,  (14) 

Mild  ~  Mild  =  diltj(£  ~  £ii)  ,  (15) 

0  —  Mifc  d  =  dik  ,j  (£  —  £i& )  , 
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(16) 


and 


(17) 


A4 it,ji  £  1,  •  .  -  ,  fc. 

Summing  (14)  -  (16)  and  using  (17)  results  in 

k 

E^d(£-&j3  =  °-  (18) 

t=0 

Similarly,  for  a  vertical  edge  segment  we  have  that 


i 

E^O?  -  Vh)3 


t=o 


=  0. 


(19) 


We  call  (18)  and  (19)  edge  conformality  conditions. 


Lemma  6.1.  If  each  and  r]je  are  different ,  then  the  dimension  of  the  solution  spaces 
corresponding  to  (18)  and  (19)  are  k  —  3  and  l  —  3,  respectively. 


Proof.  We  only  prove  the  dimension  of  the  solution  space  corresponding  to  (18).  The 
proof  for  (19)  is  similar.  Denote  dj  =  [di0  j, . . . ,  ditij, . . . ,  dikij]T  and 


'  1 

1  . 

.  1 

1 

m,  = 

& 

£ik 

£*fc+l 

’*2 

■  f2 

s*fc+i 

L 

63  • 

•  £3 

£3 

s*fc+i 

Then  (18)  is  equivalent  to  the  linear  system 


(20) 


m/Cl,  =  0. 


(21) 


The  rank  of  matrix  nx,  is  4  since  k  >  4  for  every  horizontal  edge  segment  in  Text(S). 
The  assertion  follows  from  the  rank-nullity  theorem.  □ 


Assembling  all  horizontal  and  vertical  edge  conformality  conditions  into  a  global  sys¬ 
tem  generates  the  global  conformality  condition  for  Sext-  In  other  words, 


MD  =  0  (22) 

where  DT  =  [di,d2,  ■  ■  ■ ,  dnex t]T  is  a  column  vector  of  all  vertex  cofactors  in  Text{5)  and 
M  is  a  4 nseg  x  next  real  matrix.  Each  edge  conformality  condition  corresponds  to  a 
submatrix  consisting  of  4  rows  of  M  and  each  vertex  cofactor  corresponds  to  a  column 

of  M. 

Lemma  6.2.  The  dimension  of  Sext  is  the  nullity  of  M,  i.e.,  the  dimension  is  next 
minus  the  rank  of  M. 

Proof.  Since  the  continuity  constraints  in  Sext  have  been  converted  into  the  linear  system 
in  (22),  the  dimension  of  Sext  is  the  dimension  of  the  null  space  of  M,  i.e.,  the  dimension 
is  next  minus  the  rank  of  M.  □ 
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6.3.  Simplifying  the  constraint  matrix,  M,  and  Tea;t((5) 


It  is  possible  to  simplify  the  constraint  matrix,  M,  and  the  topology  of  the  extended 
T-mesh,  Text(5)  such  that  the  null  space  of  M  is  undisturbed.  To  remove  a  vertex 
from  Text(S)  means  we  delete  the  corresponding  column  from  M  and  to  remove  an  edge 
segment  from  Tea;t(<5)  means  we  delete  the  appropriate  submatrix  from  M.  We  form 
the  reduced  constraint  matrix  M  by  removing  the  eight  edge  segments  e\,  elf,  ehh 

nseg  1 

ehh  ,  and  eY,  ef,  elv  ,  ,  and  sixteen  corner  vertices  V-2  -2,  u_i  -2,  V-2  -1,  V-i  -1 

nseg  seg  x  ,Lseg 

and  V-2,d+2,  V-l,d+2:  V-2,d+3,  V-i^d+3,  an(l  Vc+2,d+2,  W+3,d+2 ,  Vc+2,d+3,  Vc+3,d+3i  and 
uc+2,_2,  ^c+3,-2,  Vc+2,-1,  Vc+3,-i  from  Text(6).  We  denote  the  T-mesh  after  the  removals 
by  T ext{d)  and  the  number  of  vertices  and  edge  segments  in  T ext(d)  by  next  and  nseg, 
respectively  Figure  7  shows  the  simplified  extended  T-mesh  Text(i5)  corresponding  to 
the  extended  T-mesh  in  Figure  4b.  The  vertices  and  edge  segments  which  remain  after 
the  removal  process  have  corresponding  entries  in  M. 


Figure  7:  The  simplified  extended  T-mesh  Text(<5)  corresponding  to  the  extended  T-mesh  in  Figure  4b. 


Lemma  6.3.  The  dimension  of  the  null  space  of  M  is  the  same  as  that  for  M. 

Proof.  The  vertex  cofactors  which  correspond  to  the  removed  corner  vertices  can  be 
uniquely  determined  by  applying  (18)  to  the  four  horizontal  removed  edge  segments  or  by 
applying  (19)  to  the  four  vertical  removed  edge  segments.  To  establish  the  result  we  need 
to  show  that  the  constraints  corresponding  to  the  four  vertical  removed  edge  segments 
can  be  derived  from  the  constraints  corresponding  to  the  four  horizontal  removed  edge 
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segments.  We  have  that 


°=£ 


m= 0  L^— 0 


(v-Vjt) 


l 


=  £ 

^=0 


/e 


^  1  dim,je(£  £im) 


3  (£-&j3 

3  {V-Vje)3 


.171= 0 


=  £ 


£=0,1, 2-1, 2  Lm=0 


£dw<(e-&J3 


(v-riuf- 


(23) 

(24) 

(25) 


Equation  (23)  involves  the  sum  of  all  edge  conformality  conditions  for  the  horizonal  edge 
segments.  Equation  (25)  holds  because  the  linear  systems  for  the  other  vertical  edge 
segments  are  satisfied.  Since  ( p  —  Pje),  £  =  0,1,1  —  1,1,  form  a  basis  for  a  linear  space  of 
polynomials  with  degree  less  than  four,  X)m=o  <i;m ,je  (£  —  £im)3  =  0,  for  m  =  0, . . . ,  k.  In 
other  words,  the  constraints  for  the  four  vertical  removed  edges  segments  can  be  derived 
from  the  other  constraints.  □ 

Lemma  6.4.  If  M  has  full  column  rank,  then  the  dimension  of  Sext  is 

dim  Sext  =  na  +  n+  +  n~  (26) 

where  na  is  the  number  of  active  vertices  in  T (6)  and  n+  and  n~  are  the  number  of 
crossing  and  overlap  vertices,  respectively,  in  Text  (6) . 

Proof.  Since  there  are  next  and  nseg  vertices  and  edge  segments,  respectively,  in  Text(6), 
M  is  a  4 nseg  x  next  matrix.  And  M  has  full  column  rank  the  dimension  of  Sext  is 
next  ~  4 nseg.  As  every  edge  segment  in  Tea;t(5)  has  exactly  four  extended  vertices  (see 
Section  3.2)  and  these  four  extended  vertices  are  not  extended  vertices  for  any  other  edge 
segment,  the  number  of  extended  vertices  in  T ext(6)  is  4mseg.  Thus, 

dim  Sext  =  next  -  4mseg  =  na  +  n+  +  n~ .  (27) 


□ 


7.  The  dimension  of  l~ext  and  Teiem 

We  now  prove  a  fundamental  dimension  result  which  is  written  in  terms  of  the  topol¬ 
ogy  of  a  T-mesh.  This  result  is  fundamental  in  establishing  the  nesting  results  in  Sec¬ 
tion  8. 

Theorem  7.1.  For  any  T-mesh  T ,  the  dimension  of  Sext  is 

dim  Sext  =  na  +  n+  +  n~ .  (28) 

Proof.  According  to  Lemma  (6.4),  we  can  establish  the  result  if  we  can  show  that  M  has 
full  column  rank.  Since  every  extended  vertex  in  Tea;t(<5)  is  an  extended  vertex  in  exactly 
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one  edge  segment,  the  matrix  M  has  more  columns  than  rows,  i.e. ,  next  >  Anseg.  An 
appropriate  partition  of  the  linear  system  of  constraints,  M  D  =  0,  is 


[  Mi 


=  0 


(29) 


where  Mi  is  a  4nseg  x  4nseg  matrix  and  Mi  is  a  4 nseg  x  ( next  —  4 nseg)  matrix,  Di  = 
[d\, ...,  d^seg]T  is  a  vector  of  the  first  4hseg  vertex  cofactors,  and  D2  =  [dj, . . . ,  ^ }T 

is  a  vector  of  the  remaining  vertex  cofactors.  Thus,  the  problem  reduces  to  finding  an 
appropriate  ordering  of  edge  conformality  conditions  and  vertex  cofactors  such  that  Mi 
is  full  rank.  There  are  two  cases  to  consider: 


Case  1:  Text(i5)  has  no  circles 

Since  Text(5)  has  no  circles  it  is  possible  to  create  a  partial  ordering  of  the  edge 
segments  such  that  comes  before  e,j  if  no  boundary  vertex  in  e,  is  an  interior 
vertex  in  edge  segment  eg.  Then,  using  this  ordering,  for  *  =  1,  2, ... ,  hseg,  we  place 
the  four  extended  vertex  cofactors  and  edge  conformality  conditions  corresponding 
to  edge  segment  e,  in  rows  4(«— 1)  +  1  through  4(«  — 1)+4  of  Di  and  M,  respectively. 
The  interior  vertices  of  edge  segment  ti  can  be  placed  anywhere  in  D2.  As  a  result, 
the  matrix  Mi  is  in  upper  block  triangular  form  and  according  to  Lemma  6.1  each 
diagonal  block  4x4  matrix  is  full  rank,  thus  matrix  Mi  is  obviously  of  full  rank. 


Case  2:  Text(i5)  has  circles 

If  Text{5)  has  circles  the  matrix  M  is  more  complex.  For  a  circle,  ei,i  =  1, . . .  ,c 
with  circle  vertices  Vi ,  the  edge  conformality  conditions  and  vertex  cofactors  as¬ 
sociated  with  the  circle  can  be  arranged  to  create  a  circulant  block  matrix  if  we 
place  the  edge  conformality  conditions  corresponding  to  the  four  extended  vertex 
cofactors  for  in  rows  4 (*  —  1)  +  1  through  4 (*  —  1)  +  4  with  the  edge  conformal¬ 
ity  condition  corresponding  to  the  vertex  cofactor  for  the  circle  vertex  in  row 
4 {i  —  1)  +  1.  The  circulant  block  matrix  has  the  following  form, 


mi  ni  0  ...  0 

0  m2  112  •••  0 


(30) 


0  mc_i  nc_i 

nc  0  ...  0  mc 

where  m,,  i  =  1,.  ,.,c,  is  a  4  x  4  matrix  corresponding  to  the  edge  conformality 
conditions  for  the  four  extended  vertices  in  circle  edge  segment  and  ,  i  = 
1, . . . ,  c  —  1,  is  a4x4  matrix  corresponding  to  the  edge  conformality  conditions  for 
circle  edge  segment  et  where  only  the  first  column  is  non-zero. 


Through  row  reduction  we  can  transform  the  circulant  block  matrix  into  an  almost 
triangular  matrix 

I  ni  0  ...  0  ' 

0  I  fl2  0 


1 
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m, 


0 

0 


0 

0 


I 

0 


(31) 


where  I  is  a  4  x  4  identity  matrix,  n  i,  i  =  1, . . . ,  c  —  Iisa4x4  matrix  corresponding 
to  the  edge  conformality  conditions  for  circle  edge  segment  e*  where  only  the  first 
column  is  non-zero,  and  mc  is  a  4  x  4  matrix  where  the  first  row  and  column  may 
be  zero.  If  so,  we  call  the  circle  vertex,  vc,  a  vanished  vertex  and  the  circle  edge 
segment,  ec,  a  vanished  edge  segment.  We  can  now  make  the  following  observations 
about  circles: 

1.  A  circle  has  at  most  one  vanished  vertex  which  depends  on  the  ordering  of 
the  edge  conformality  conditions.  In  other  words,  if  we  arrange  the  circle  edge 
segments  as  et+i,  . . . ,  en,  e±,  . . . ,  e,,  then  the  vanished  vertex  will  be  Vi,  and 
the  vanished  edge  segment  is  e,;. 

2.  If  T ext{5)  has  hc  circles  then  there  are  at  most  hc  vanished  circle  vertices. 
Since  each  circle  has  at  least  two  horizonal  edge  segments  and  each  edge 
segment  belongs  to  at  most  four  circles  we  have  that  2 hgeg  >  nc,  where 

is  the  number  of  horizontal  edge  segments  in  Tea,t(<5).  Thus,  we  can  create  a 
selection  set  of  horizonal  edge  segments,  e±,  . . . ,  e^,  nt  <  hc  such  that 

each  edge  segment  in  the  set  is  a  vanished  edge  segment,  each  belongs  to 
different  circles,  and  each  has  at  most  two  vanished  extended  vertices. 

For  the  edge  segments  which  are  not  part  of  the  selection  set  we  can  create  a 
partial  ordering  as  described  in  Case  1.  We  place  the  vertex  cofactors  and  edge 
conformality  conditions,  corresponding  to  each  edge  segment  in  the  partial  order¬ 
ing,  into  D  and  M  as  described  in  Case  1.  Then,  for  each  edge  segment  in  the 
selection  set,  we  first  place  the  vertex  cofactors  corresponding  to  two  active  vertices 
into  Di,  followed  by  two  non- vanished  extended  vertices.  The  corresponding  edge 
conformality  conditions  are  placed  in  Mi.  By  construction  Mi  has  the  following 
form 


’  m™t+l 

0 

0 

0 

* 

rxin 

ii'seg 

* 

Cl 

0 

0 

c 

0 

0 

0 

0 

cnt 

The  4x4  matrix,  mi,  i  —  fit  +  1, . . .  ,hseg,  corresponds  to  the  edge  conformality 
conditions  for  the  four  boundary  vertices  in  the  (non-circle)  edge  segment  e*.  The 
4x4  matrix,  c,;,  i  =  1  ,...,ht,  corresponds  to  the  edge  conformality  conditions 
for  the  two  non-vanished  extended  vertices  and  two  active  vertices  in  circle  edge 
segment  ej.  The  4 ht  x  4 (hseg  —ht  + 1)  matrix,  C,  corresponds  to  the  vertex  cofactor 
coupling  in  the  circles.  The  *  means  that  the  corresponding  part  of  the  matrix 
may  be  non-zero  but  is  not  important  for  our  developments.  We  now  perform  row 
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reductions  to  reveal  the  rank  of  Mi.  The  first  pass  of  row  reductions  yields 


"  mnt  +  l 

0 

0 

0 

* 

* 

I 

0 

0 

c 

0 

0 

0 

0 

I 

(33) 


A  second  pass  of  row  reductions  eliminates  matrix  C  and  produces  the  following 
upper  block  triangular  matrix 


'  m^+1  =1= 

0 

* 

0  0  m„ 

ILseg 

Ci  0  0 

0 

0  0 

0  0  c  nt 

(34) 


The  4x4  matrices,  c*,  i  =  1, . . .  ,fit,  are  diagonal  matrices.  Notice  that  only  the 
information  of  circle  Ci  is  required  to  eliminate  the  edge  conformality  conditions  for 
edge  segment  e,.  Since  no  vertex  cofactors  corresponding  to  vanished  vertices  are 
in  Di  no  entries  of  c i  are  zero.  Thus,  Mi  is  full  rank. 

We  can  now  conclude  that  M  is  full  column  rank.  □ 

Corollary  7.2.  For  any  T-mesh  T,  the  dimension  of  Text  is 

dvcctText  =  na  +  n+  +  n~ .  (35) 

Proof.  Apply  Theorem  5.1.  □ 

Corollary  7.3.  Given  T1  and  T2  such  that  T\xt  C  T2xt,  then 

dim  7 eL  =  dim  T^xt  +  ndiff  =  na  +  n+  +  n~  +  ndiff  (36) 

where  ndl^  is  the  number  of  vertices  in  T2xt  which  are  not  in  T \xt . 

Proof.  Since  the  number  of  edge  segments  in  both  T-meshes  is  the  same  we  can  place 
the  vertex  cofactors  corresponding  to  the  vertices  in  T2xt  (5)  which  are  not  in  T \xt  {8)  into 
D2.  This  implies  that  Mi  is  identical  for  both  T-meshes.  This  implies  that  dimS2^  = 
dimiSg^t  +  ndlff  =  na  +  n+  +  n~  +  ndl^ .  Then  apply  Theorem  5.1  to  arrive  at  the 
result.  □ 
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8.  The  nesting  behavior  of  T-spline  spaces 

We  can  use  the  dimension  results  derived  for  Sexti  Seiem,  Text ,  and  Teiem  to  prove 
several  fundamental  nesting  properties  for  T-splines  spaces.  Specifically,  we  prove  Theo¬ 
rem  8.9  and  Corollary  8.10  which  serve  as  the  theoretical  foundation  for  analysis-suitable 
local  refinement  [3] . 

Lemma  8.1.  For  any  T-mesh ,  T.  Text  Q  Tiem 

Proof.  This  immediately  follows  from  the  fact  that  Text  C  T eiem.  □ 

Lemma  8.2.  For  any  T-mesh ,  T,  T(8)  C  Teiem- 

Proof.  This  immediately  follows  from  the  fact  that  every  blending  function  N  a  £  T(8) 
is  in  Tei  em  •  □ 

Lemma  8.3.  There  exists  T-meshes  such  that  T(8)  ^  Text- 

Proof.  Consider  the  T-mesh  and  superimposed  extended  and  elemental  T-meshes  in  Fig¬ 
ure  8.  The  extended  T-mesh,  T ext(5),  is  denoted  by  the  dashed  lines  while  the  elemental 
T-mesh,  Teiem(8),  is  denoted  by  the  bold  shaded  lines.  In  this  case,  T{S)  C  Teiem  and 
T(8)  £  Text-  □ 


Figure  8:  A  T-mesh  with  the  extended  and  elemental  T-meshes  superimposed  on  one  another.  In  this 
case  T(S)  C  Teiem  and  T($)  £  Text •  The  extended  T-mesh  Text(<^)  is  denoted  by  the  dashed  lines  while 
the  elemental  T-mesh  Te^em(<5)  is  denoted  by  the  bold  shaded  lines. 


Definition  8.4.  A  polynomial  T-mesh  has  an  extended  T-mesh ,  T ext,  with  the  following 
properties: 

•  There  are  no  crossing  vertices. 

•  extir)  —  Teiem(e). 
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A  T-spline  space  with  a  polynomial  T-mesh  is  called  a  polynomial  T-spline  space 


To  illustrate  the  definition  consider  the  T-meshes  in  Figure  9.  The  T-meshes  on  the 
right  are  e-offset  T-meshes  corresponding  to  the  T-meshes  on  the  left.  In  all  cases,  the 
extended  and  elemental  T-meshes  are  superimposed  on  one  another  with  the  dashed  lines 
denoting  an  extended  T-mesh  and  the  bold  shaded  lines  denoting  an  elemental  T-mesh. 
The  T-mesh  in  Figure  9a,  on  the  left,  is  a  polynomial  T-mesh  since  Text(e)  =  Te;em(e) 
as  shown  on  the  right.  The  T-mesh  in  Figure  9b,  on  the  left,  is  not  a  polynomial  T-mesh 
since  Text(e)  ^  Te;em(e)  as  shown  on  the  right.  Notice  that  this  subtle  but  important 
difference  is  not  discernible  when  the  T-meshes  are  not  e-offset. 


(a)  A  polynomial  T-mesh 


(b)  Not  a  polynomial  T-mesh 


Figure  9:  The  T-mesh  on  the  top  left  is  a  polynomial  T-mesh  while  the  T-mesh  on  the  bottom  left  is 
not  a  polynomial  T-mesh.  This  can  be  seen  by  examining  the  corresponding  e-offset  T-meshes  on  the 
right. 


Definition  8.5.  An  analysis- suitable  T-mesh  has  an  extended  T-mesh,  T ext,  which  does 
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not  have  any  edge  or  face  extension  intersections.  A  T-spline  space  with  an  analysis- 
suitable  T-mesh  is  called  an  analysis-suitable  T-spline  space. 

Lemma  8.6.  An  analysis- suitable  T-spline  space  is  a  polynomial  T-spline  space. 

Proof.  By  definition,  Text,  corresponding  to  an  analysis-suitable  T-spline  space,  has  no 
face  extension  intersections.  To  prove  that  Text(e)  =  Te;em(e)  see  [17],  Lemma  9.  □ 

Lemma  8.7.  For  a  polynomial  T-mesh,  dim  7^4  =  na  +  n~ . 

Proof.  This  follows  from  Corollary  7.2  and  the  fact  that  n~  =  0  for  a  polynomial  T- 
mesh.  □ 

Lemma  8.8.  For  a  polynomial  T-mesh,  T{8)  C  Text- 

Proof.  According  to  Lemma  8.2,  T{8)  C  Teiem-  If  T  is  a  polynomial  T-mesh,  Text  = 
Telem  -  I— 1 

Theorem  8.9.  Given  two  linearly  independent  polynomial  T-meshes,  T1  and  T2 ,  if 
Tlxt  Q  Tlx t,  then  T1  CT2. 

Proof.  It  is  obvious  that  Tf}xt  C  Text-  If  both  T-meshes  have  no  zero  knot  intervals  we 
consider  the  following  three  cases: 

1.  T \xt  and  T2xt  have  no  overlap  vertices. 

In  this  case,  T1  =  TLt  an(I  "T2  =  "Qxti  so  T1  C  T2  ■ 

2.  T lxt  has  overlap  vertices  and  T2xt  does  not. 

In  this  case,  T1  C  Tfxt  and  T2  =  T^xt,  so  T1  QT2. 

3.  T2ext  has  overlap  vertices 

Form  Tx{e)  and  7^(e).  According  to  Case  1,  Tx{e)  C  T 2(e).  Now  let  e  — >  0. 
Otherwise,  since  Tlext{8)  C  T2xt{8)  we  have  that  Tl{5)  C  T2(<5).  Now  let  <5  — >-  0.  □ 

Corollary  8.10.  Given  two  analysis-suitable  T-meshes,  T1  and  T2,  if  T\xt  C  T2xt,  then 

T1  c  T2. 

Proof.  This  follows  from  Theorem  8.9  and  the  fact  that  analysis-suitable  T-meshes  are 
linearly  independent  (see  [16],  Theorem  18).  □ 

9.  Conclusion 

We  have  established  the  nesting  behavior  of  analysis-suitable  T-spline  spaces.  This 
provides  a  theoretical  foundation  for  the  efficient  and  highly  localized  local  refinement 
algorithm  presented  in  [3]  and  its  use  in  isogeometric  analysis.  Additionally,  we  have 
derived  a  dimension  formula  for  smooth  polynomial  spline  spaces  defined  over  the  Bezier 
mesh  of  a  T-spline  and  demonstrated  its  utility  in  developing  the  theory  of  T-spline 
spaces. 
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